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Abstract: We use the AdS/CFT correspondence to determine the out-of-equilibrium 
production rate of dileptons at rest in strongly coupled N = 4 Super Yang-Mills plasma. 
Thermalization is achieved via the gravitational collapse of a thin shell of matter in AdSs 
space and the subsequent formation of a black hole, which we describe in a quasistatic ap- 
proximation. Prior to thermalization, the dilepton spectral function is observed to oscillate 
as a function of frequency, but the amplitude of the oscillations decreases when thermal 
equilibrium is approached. At the same time, we follow the flow of the quasinormal spec- 
trum of the corresponding U(l) vector field towards its equilibrium limit. 
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1 Introduction 

Experiments conducted at RHIC [1] and at the LHC [2] suggest that the matter produced 
in heavy ion collisions can be identified as strongly coupled quark-gluon plasma, behaving 
essentially like a strongly coupled liquid. The AdS/CFT correspondence [3-5], on the 
other hand, is a powerful tool for studying the properties of strongly interacting, large- N 
field theories at finite temperature and density; for recent reviews see e.g. [6, 7]. In the 
AdS/CFT setting, thermalization on the field theory side is related to the classical gravity 
description of black hole formation in a 5-dimensional spacetime. 

While holographic thermalization is in general a very complicated problem (discussed 
e.g. in [8-10]), an interesting approach was taken in [11, 12], where the authors studied the 
quasistatic formation of a black hole through the gravitational collapse of a thin shell of 
matter in AdS space, falling in the radial direction. Amongst other things, they determined 
the retarded Green's function of a massive scalar field living on the boundary of AdS space, 
as the shell radius decreased towards the Schwarzschild radius and the field theory system 
approached equilibrium. In [13, 14], this picture was further extended by analyzing the 
correlators of metric perturbations. 

In the paper at hand, we consider the 72-current correlator in strongly coupled M = 4 
supersymmetric Yang-Mills theory (SYM) in the large N limit, having in mind the spectral 
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function related to the production of dileptons at rest [15]. On the gravity side, this 
amounts to introducing an effective U(l) gauge field in AdSs space and determining the 
corresponding spectral function as well as the quasinormal spectrum. In contrast to earlier 
works on the subject, we want to perform the calculation in an out-of-equilibrium setting, 
and to this end follow references [11, 12] by introducing to the system a thin, spherical shell 
of matter that undergoes gravitational collapse in the fifth (radial) dimension of the AdS 
space. When the radius r s of the shell approaches the Schwarzschild radius dictated by 
its mass, the dual field theory system undergoes thermalization. In the present paper, we 
restrict ourselves to a limit where the motion of the shell can be considered quasistatic, 
and hence the production rate of dileptons is available using methods discussed e.g. in [16]. 
The (QCD version of the) quantity we compute is measurable in real world experiments 
and is moreover a very interesting probe of the initial stages of a heavy ion collision; for 
reviews, see e.g. [17-19] and references therein. 

Our paper is organized as follows: In section 2, we introduce the setup in AdS5 space 
by writing down the corresponding metric and discussing the matching conditions for dif- 
ferent types of fields at the collapsing shell. In Section 3, we then determine the relevant 
holographic Green's functions and relate them to dilepton production, while in section 4 we 
draw our conclusions. In appendix A, we finally discuss the case of a massive scalar field in 
3-dimensional AdS3 space, and in appendix B the validity of the Fluctuation Dissipation 
Theorem (FDT) in our setup. 

After the first version of our paper was posted on the arxiv, several interesting and 
closely related papers on holographic thermalization appeared. In [20], the authors analyze 
the falling shell setup and manage to go beyond the quasistatic approximation in their study 
of unequal-time correlators. Ref. [21] on the other hand studies holographic thermalization 
in the presence of a finite chemical potential, while ref. [22] attempts to derive a holographic 
FDT valid out of thermal equilibrium. 



2 Setup and matching conditions 

Our aim is to follow the idea of Danielson et al. [11, 12] and describe the thermalization 
process in strongly coupled, large- N N = 4 SYM plasma via the gravitational collapse of an 
infinitely thin shell of matter in a background, which is asymptotically AdSs. Analogously 
to Birkhoff's theorem, we know that outside the shell the metric is given by a black hole 
solution, i.e. the Schwarzschild solution in AdSs; similarly, the metric inside the shell is 
given by the flat AdSs metric. Thus, we can write the metric in our setup in the form 

ds 2 = -f{ r )dt 2 + ^+r 2 d* 2 , (2.1) 
fir) 

where x stands for the coordinates of Euclidean 3-space and we have set the curvature 
radius of AdS space to unity. If the shell resides at r = r s , then the function /(r) is given 
by 

r( x ff-(r) = l + r 2 , for r < r s 

f( r > = \ f , n -i m 2 , 2 t ^ ' ( 2 - 2 ) 
/+(r) = 1 - + r z , for r > r s 
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where the parameter m is related to the mass of the shell. To match the two (in principle 
unrelated) coordinate patches, we choose the obvious condition x|+ = x|_. From this, it 
follows that the r-coordinate is continuous over the boundary, while the time coordinate t 
has a discontinuity at the shell. Finally, in the limit where the mass of the shell is large 
enough so that the corresponding horizon radius 3> 1, we can relate the parameter m 
to r/j and the Hawking temperature of the black hole T through 

m = r \ = ir 2 T 2 . (2.3) 

In the presence of a black hole, i.e. once the thermalization process has ended, the parameter 
T is identified with the equilibrium temperature of the field theory; in the rest of this paper, 
we will, however, refer to this definition of T even when the shell resides at a radius r s > r^. 

We are interested in finding matching conditions for various types of fields at r = r s , 
and to this end, we need the normal vector of the shell. Let the coordinates of the shell be 

t = t a (r) , r = r s (r) , 

where the parameter r is the proper time of the shell. From the condition u^u^ = —1, we 
find 

-/ ± (r,« + ^ = -l, (2.4, 

where we have denoted a derivative with respect to r by a dot; from here, it also follows 
that 

; y/f±(r 8 ) + r 2 s 

ts± = Mrs) ■ (2 - 5) 
Next, we introduce the normal vector of the shell, n^, and require it to be normal to the 
four-velocity of the shell as well as properly normalized. This produces 



n t i s + n r r s = , (2.6) 
-j^- ) +f±(r s )n 2 r = l, (2.7) 

which after some algebra leads us to the simple result 

[n„] = (-r s , 4,0,0,0) . (2.8) 

Note that while r s is the same on both sides of the shell, t s is not, but has a discontinuity 
given by eq. (2.5). 

To determine the trajectory of the shell, r s (r), one should next specify its energy con- 
tent, after which one would proceed to solve its geodesic equation. While this information 
is in principle crucial to disentangle the time evolution of the thermalization process, we 
will in the present work restrict ourselves to the quasistatic limit, where the frequency of 
the modes we study is assumed to be much larger than the inverse time scale related to the 
falling of the shell. This produces a crucial technical simplification, and corresponds to a 
particular choice of initial conditions for the thermalization process. In section 4, we will 
inspect the validity of this approximation, and in addition relate the value of the parameter 
m to the characteristic time scale of the thermalization process. 
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2.1 Junction condition for a scalar field 



Let us first study the behavior of a scalar field (p at the shell radius, both for simplicity 
and for the eventual use of the results in appendix A. Its junction condition states that the 
covariant derivative normal to the shell should be continuous across the shell [23], i.e. 



Writing this out explicitly in coordinate space gives 



f-(r.)M. 



Or 



f+(r s )dt, 



+ yjf + (r 8 )+r1 



Or 



(2.9) 



, (2.10) 



where t± reminds us of the fact that the time coordinates are different on the two sides of 
the shell, and thus need to be evaluated separately. 

The quasistatic approximation is equivalent to assuming that the timescale related 
to the motion of the shell is longer than any other timescale of interest. We may thus 
approximate r s — > 0, which gives us 



Or 



dr 



Assuming further the continuity of (f)(t) across the shell, (ft-(t- 
relation 



dt. 
dL 




I 



1 



(2.11) 

<fi + (t + ), and using the 

. UJj_t_ 

(2.12) 



dt-e vT^l 



this leads to the identification w_ 
conditions 



u+/\fjm as well as to the Fourier space matching 



Mw-) = Vfrn<t>+(u+), (2-13) 

= f m <t>' + (2.14) 

where the prime denotes a derivative with respect to r. 

2.2 Junction condition for a vector field 

To find the junction condition in the vector field case, we loosely follow the derivation given 
in [23]. We first rewrite the metric in the form 

ds 2 = _ fdt 2 + j! +r 2 dx 2 5 (215) 



f± 



where the /-functions are given by 

/ = 

/ = 





for 


r < r s 




for 


r > r s 




for 


r <r s 








<\ 


for 


r > r s 



(2.16) 
(2.17) 
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This is achieved by first approximating r> 1, and then rescaling time in the r < r s patch. 

In the above coordinate system, the equation of motion for an electric field in the z 
direction, E = iF tz , is given by (see also eq. (3.3) of the following section) 




(2.18) 



Introducing the tortoise coordinate 



dr* 
dr 




1 



(2.19) 



this simplifies to the form 



d 2 r E + u?r 2 E = 0. 



(2.20) 



From this equation, we get 



d r ,E\± 




uj 2 r 2 E dr. 



(2.21) 



where the right hand side furthermore vanishes in the limit of an infinitesimally thin shell 
(e — > 0). Going finally back to our original coordinate system and taking into account the 
jump in frequencies, we obtain from here a result very similar to the scalar field one, 



3 Holographic dilepton production 

In this section, we want to determine the production rate of dileptons at rest in the holo- 
graphic setup introduced in section 2. We begin this in section 3.1 by discussing, how the 
dilepton production rate is related to the retarded Green's function of a photon field, and 
then proceed to write down the equations of motion for a U(l) vector field in the presence 
of a falling shell in section 3.2. After this, we evaluate the necessary retarded correlators 
in section 3.3, while section 3.4 is devoted to a study of the quasinormal spectrum of the 
gauge field, providing another way to follow the thermalization process. In section 3.5, we 
finally analyze the large \u>\ limit of the correlator in more detail, applying to it the WKB 
approximation. 

3.1 Dilepton production rate from Green's functions 

It is a textbook exercise to show that the differential production rate of dileptons with 
four- momentum Q is directly proportional to the photon Wightman function H<(Q) (see 






(2.23) 



e.g. [24]) 



d 4 Q 



dr 



mf n<„(Q) 

24tt 4 Q 2 
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where a is the fine structure constant. Just as in thermal equilibrium, one can show that 
in the quasistatic approximation the Wightman function can be further related to the 
corresponding spectral density Xi i- e - the imaginary part of the retarded Green's function 
Hpu- 1 Setting the external three-momentum q to zero, which implies that the produced 
dileptons are at rest and that the longitudinal and transverse parts of U pu agree, 2 we obtain 

rrn<» = -2n(w)ImIl£( W ) = ^(u), (3.2) 

where n{oj) stands for the usual Bose-Einstein distribution function. Our task thus becomes 
that of evaluating the retarded Green's function of a U(l) gauge field in the gravity picture. 

3.2 Equations of motion and analytic solutions 

The equation of motion of a U(l) vector field in curved spacetime reads 

' dp [V^99 w 9 ua F pa ] = , F pa = d p A a - 8 a A p , (3.3) 



V~9 

where the metric is in our case given by eq. (2.1). Once again, we approximate r>l and 
write the mass in terms of the horizon radius r^. To simplify the otherwise cumbersome 
expressions, we also introduce the commonly used variables 

u = -| and z = 1 — u , 

which both vary between and 1. In the following, we will repeatedly interchange between 
the three radial variables; it should, however, always be clear from the context, which one 
is being used at a given time. 

The gauge invariant component of the field strength tensor we are interested in can be 
chosen as Ft z (t, tt,x). Given in terms of the Fourier components of the vector potential, 
Apty, u,x) = Jq e~ %U3t+icl '' x -A ti (cj, u, q), this quantity reads [15, 16] 

Ff Z (u, u, q = 0) = — iu)A z (u, u) = —iE(to,u), (3.4) 

where E stands for the electric field in the z direction and we have set the three-momentum 
q to zero. A short calculation gives for the equation of motion of E 

d\E + ^ld z E + {1 f z)f2 E = , (3.5) 

where we have introduced the dimensionless variable 

to to 

to = 



2r h 2ttT 



1 This relation, a special case of the Fluctuation Dissipation Theorem, can be derived by generalizing 
the arguments of Ref. [25] to our present case, where instead of a black hole we have a static shell sitting 
at a radius r s > r^. For further details, see appendix B. 

2 This in particular means that it suffices to consider a single Green's function II from this point onwards. 
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This equation can be identified with the Riemann differential equation [26], 



diE + 



1 1 

+ 



d z E + 



co 2 /2 u 2 /2 



E 



z ' z — 2\ " z z — 2 

which has two linearly independent solutions, 

iui (l+i)£> a ( 1 + % 

E iQ (u,z) = z~r(l - z)^(2 - z)-2 2 F 1 M - — u, 
E out (co, z) = zf(\- z) (2 ^(2 - z)-f 2 Fi (1 - ^-u, 



z(*-l)(«-2) 



0, 



1 + i 



w; 1 — id); 



cj; 1 + iuj] 



2(z-i; 

z 

2(^-1) 



(3.6) 

(3.7) 
(3.8) 



that can be expressed in terms of the hypergeometric functions 2 -Pi (see also [16]). Here, 
i?i n is the so-called infalling solution, obeying the incoming wave boundary condition at 
the horizon, z — > 0, while E ont satisfies the outgoing boundary condition in this limit. In 
the limit z — > 0, or r — > r^, the two solutions can be approximated by the expressions 



E in (uj,z^ 0) 



out 



(3.9) 



In the presence of a black hole, the physical solution is the one obeying the infalling 
boundary condition, as classically nothing can escape from inside the event horizon [27]. 
In our case, waves can, however, be reflected from the shell and escape to the boundary, 
so outside the shell we must write the general solution as a linear combination of the two 
independent solutions (here and in the following, we use u to denote the frequency outside 
the shell), 

-Eoutside^, r) = c + E in (u, r) + c-E out (uj, r) . (3.10) 

Inside the shell, r < r s , the solution is on the other hand obtained from the / — > 1 limit of 
eq. (3.5), giving 



E 



inside 



(w,r) 



1 



J1 



^ ^inside ^ 



+ iY l 



^ ^inside ^ 



r-+Q 1 



^inside A 



f V ^inside 



(3.11) 



where again upside = ^/VJm, with f m defined in eq. (2.12). The matching conditions of 
eq. (2.22) then lead to the result 



c_(r s ) E in (u), r)d r E insidc (uj, r) - \ZjmE imjide (uj , r)d r E in (u, r) 



(3.12) 



c+(r s ) E out (u, r)d r E- msidc (uj, r) - s/jmE inside (u> , r)d r E out (uj, r 
which we will use frequently below. For real values of u, this ratio satisfies the relation 

ct(r s ) 



c-(r. 



(3.13) 



that will prove rather useful in the following. 
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3.3 The retarded Green's function and spectral density 

The retarded correlator (Green's function) of an operator living on the boundary of AdS 
space can be determined using the standard AdS/CFT prescription, developed and thor- 
oughly explained in [27, 28]. The wave function of the corresponding bulk field is first 
written in terms of a power series expansion around the singular point (s) of its equation 
of motion (in our case -^outside at u = 0), 

Outside (w, u) U =° A(u)u A - [1 + . . .] + B(co)u A + [1 + ...}, (3.14) 

where A± are the related characteristic exponents. For -Eoutsidei we have A+ = 1 and 
A_ = 0, and hence the solution takes the form 

E outside (oo,u) A(u) [1 + h(cj)uln(u) + ...] + B(u) [« + ...] . (3.15) 

The retarded Green's function of the boundary field is then given by the expression 

which can equivalently be expressed in terms of its thermal limit times a correction function, 

1%) = n thermal (w) x H(u, r s ) , (3.17) 

where we have explicitly indicated the fact that the function H depends on the radius of 
the shell, and where obviously H — > 1 as r s — > r^. 

For our electric field, the thermal limit of the retarded Green's function can be read 
off from the literature (see e.g. [16]), 

n 2 t 2 n 

j-j-thermal^^ _ 1 - ' L > " 



8 An 
N 2 T 2 



(,3.18) 



while the function H obtains the form 



C— 




Bout 


c+ 


r a 


Bin 


C— 




Aout 




r s 





1 + 

H(u,r s ) = . (3.19) 

1+ " 

Here, ip is the logarithmic derivative of the gamma function, while 7 ~ 0.5772 is the Euler- 
Mascheroni constant. The imaginary part of the thermal Green's function is identified with 
the spectral function, which represents the thermal production rate of dileptons at rest. It 
is a smooth function of u > and reads 

x the r mal (a;) = _ 2 j m prthcrmal = ^T 2 TTU 2 sinhfru) 

4 cosh(7ra)) — cos(ttlj) 
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0.0 0.5 1.0 1.5 2.0 



Figure 1. The spectral density x(u;,r s ) as a function of & for r s = rh — 5 (thermal case; dashed, 
black curve), r s = 5.01 (blue) and r s = 5.5 (red). Note that we use units, in which the radius of 
AdS space has been set to unity. 

Next, let us study the Green's function outside of thermal equilibrium. From an 
expansion of the infalling solution (3.7), one can read off the coefficients in eq. (3.15) as 

_ T(l -iw)exp [-!±iwln2] 

An " r(i-!±^)r(i + !^) ' (3 - 21) 



, ■ ( ^-u) + 1> {-^ l y LQj ) + 1,1 2 + 2 ■ - 1 



, (3.22) 



B m = Ain(w) — Co 2 

while their outgoing counterparts are obtained through the replacements 

A)ut(^) = Ai n (iu) -)■ -iu>) , jBout(w) = Bin(i& -»■ -iCj) ■ (3.23) 

From eq. (3.12), we further see that our results fulfill an important consistency condition: 
As the shell approaches the horizon, i.e. r s — > and f m — > 0, 

c +( r ) _ .-02-HflA ( Ts 



lim = -» 1 & — - 1 -»■ °° (3.24) 

r^r sC _(r) \ r ft / 

and we recover the thermal correlator. In fig. 1, we plot the spectral density x( w ; r s) = 
— 2ImII(u;) as a function of frequency both for the thermal and non-thermal cases. 

In order to further study the thermalization process, we next follow Lin and Shuryak 
[14] and plot in fig. 2 the relative deviation of the spectral density from its thermal limit, 

for different values of r s . As the shell approaches the horizon, the oscillations visible in 
the function R increase in frequency and decrease in amplitude until they finally vanish, 
as the medium thermalizes. For purposes of clarity, we have chosen not to display here the 
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entire range of R, which approaches -1 as to — > for all values of r s . 3 Prom figs. 1 and 2, 
we further see that as expected, the approach of the system towards thermal equilibrium 
is of the "top-down" type: The energetic modes, i.e. high mass dileptons, equilibrate first. 
This feature of strong coupling thermalization was already observed e.g. in [9, 14] (for a 
discussion of thermalization in a different context, see also [30] and references therein). 
It is to be contrasted to the behavior of weakly coupled "bottom up" scenarios such as 
[31, 32], in which the soft modes thermalize before the hard ones and one would naively 
expect the amplitude of R to increase as a function of Cj. 

3.4 Quasinormal modes 

Analyzing the quasinormal spectrum of a field living in the bulk — available from the 
poles of its retarded Green's function — offers a convenient way to study the effects of the 
plasma on the corresponding field theory excitation, as well as the thermalization process 
itself. On the field theory side, the spectrum gives the dispersion relation of the mode 
in question, which in the limit r s — > r\ t should approach the thermal result. A departure 
from equilibrium is expected to lead to deviations from this spectrum, as well as to the 
appearance of altogether new poles in the correlator. At zero three-momentum, these poles 
have the generic form 

«n = M n - h n , (3.26) 

where the real part M n is identified with the mass of the resonance and the imaginary part 
T n with the corresponding decay rate. 

Recall that we can write the full retarded correlator in a factorized form 

n(w) = n therma V) x H{u, r s ), (3.27) 

3 Note that this descrease in R occurs for values of cj that are outside of the region of validity of the 
quasistatic approximation, cf. section 4. See also [29] for discussion of a similar behavior of the spectral 
density, observed in a somewhat different setting. 



- 10 - 




-1.2 



Figure 3. The poles of the full retarded correlator for r s = 5.5 (with rt = 5). 



where n thermal and H are given by eqs. (3.18) and (3.19), respectively. The thermal cor- 
relator is known to possess a quasinormal spectrum with poles located at [16] (see also 
[33, 34]) 

u> = 27rTn(±l - i) = 2r h n(±l -i) , n = 0, 1, . . . (3.28) 

which we expect to be shifted when r s ^ r^. The approach of the out-of-equilibrium 
spectrum towards the thermal one can be analyzed using eq. (3.24): We see that for all 
modes satisfying the condition Im(u;) > — r^, the ratio c_/c + approaches as r s — > rh, and 
thus the poles approach the thermal form. For Im(u) < —rh, the quasistatic approximation 
on the other hand appears to break down, and hence the thermal modes with n > will 
not be seen in the r s — > rh limit. 

In figs. 3 and 4, we show the poles of the full retarded Green's function II for r s = 5.5 
and r s = 5.01 (again with = 5), displaying only the lower complex halfplane, as the 
retarded correlator is analytic on the upper half. We observe that the poles of n thermal 
have disappeared, as they get canceled by the zeros of H. At the same time, three new 
types of poles have appeared — a structure that is intimately connected with the matching 
conditions, as they enter in the ratio of eq. (3.12). 

The tower of poles denoted by the blue circles are similar to the ones found in [12], and 
in the limit of r s — > rh one can understand their flow analytically. Using the approximation 
of eq. (3.24) together with 



A 



out 



-ie 



2i£(ln 2+vr/4) 



(3.29) 



we obtain the simple equation 



n, 

2 



rh 



le 



-uj, 



(3.30) 
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Figure 4. The poles of the full retarded correlator for r s — 5.01 (with — 5). 



where we have defined a parameter (the physical meaning of which will become clear later) 



^echo 



1 



- + In 2 - In ( ^ - 1 
2 \r h 



(3.31) 



The solutions of this equation clearly coincide with the blue circular poles shown in figs. 3 
and 4. As r s — > (see fig. 5), the flow of the lowest pole approaches u = 0, indicating a 
branch point at this value [11]. 

Next, let us look at the series of poles on the lower left quadrant of the complex plane, 
denoted by the brown diamonds in figs. 3 and 4. As is evident from the figures, these poles 
have the property that they rapidly move towards the origin as r s — > r^. Indeed, expanding 
c_ / c+ in powers of r s — rh while keeping to = uj/ \/r s /rh — 1 fixed and using the fact that 
as uj — > 0, «4out/^in — >■ 1, we obtain the simple result that the poles of the Green's function 
are given by the zeros of the Hankel function of the first kind, 



hS 1} (Co/2) = 0. 



(3.32) 



A simple numerical exercise indeed verifies this analytic result in the r s —> limit. 

Finally, we have the single pole denoted by the red square, whose flow exhibits a 
qualitative change when Imw reaches the value — r^; cf. fig. 5. For Imw < — rh, its location 
can be described by the equations Re to = and 



ri, 



2l-Imw/(2r h ) I Th 



1 



(l-Imw/rjO/2 



A 



out 



(3.33) 



where as a function of Imw is real and 0(1) in the region of interest. However, when 
r s gets exponentially close to the horizon, the pole starts to approach the origin as 



Imw 



' 2 \r h 



1/2 
1) > 



(3.34) 



- 12 - 




Figure 5. The flow of the red and lowest blue pole of the retarded correlator on the complex 
to plane, as the shell approaches the horizon at = 5. We start with the shell being located at 
r s = 5.5 and follow the flow of the poles until r s — = 1CP 15 . 



It would be tempting to relate this behavior to that of a diffusive pole, appearing at fi- 
nite external three-momentum and reflecting the diffusive relaxation of large scale charge 
density fluctuations around thermal equilibrium [28]. So far we have, however, not imple- 
mented q in our calculation, and hence postpone this study for the future. 

3.5 The WKB approximation 

So far, we have mainly studied the spectral function and quasinormal spectrum of our 
gauge field for small or moderate frequencies. In the opposite limit of very large \lo\, it 
is possible to gain analytic understanding of the solutions, as one may solve the retarded 
correlator in the WKB approximation, denoting E(u) = l/vl — u 2 Z(u) and transforming 
eq. (3.5) into a Schrbdinger type form 



d 2 u Z(u) + 



or + u 



u(l — u 2 ) 2 

The WKB solution of this equation is given by [14, 17] 



Z(u) = . 



7WKB ! 



U 



S(u) 



-1/2 



exp 



Hui I S(u)du 




(3.35) 



(3.36) 



where S(u) =- y u ( 1 _ u2)2 . 

The full solution for Z can again be written as a linear combination of infalling and 
outgoing modes, 

Z = c+Z[ n + c-Z ou t . 
Close to the horizon, i.e. at u — > 1, eq. (3.35) reduces to 



Z"(u) + 



4(1 -u) 



-Z(u) = 0. 



(3.37) 
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the solutions of which can be written as 

Z in (u) = (l-u)^)/ 2 . 

out 

Taking the u — >• 1 in eq. (3.36), we on the other hand see that 



r^WKB 



(it) \/2(l — u) 1 ^ exp 



7T 



±iu [ - +ln2 



(3.38) 



(3.39) 



implying that the two WKB solutions ^y KB can directly be identified as infalling and 
outgoing modes. This leads us to the compact result 

Z 



WKB/ \ c + r •* 1 yWKB/ \ i c ~ 1 yvi 

(u) = — ^= exp [— zwaoj (tt) H — 7= exp [iujao\ Z_ 



WKB 



V2 



V2 



(«) 



(3.40) 



where ao = vr/4 + In 2. 

Next, we look at the opposite limit of u — )■ 0, in which we wish to again relate the 
WKB solutions to the exact ones. In this limit, eq. (3.35) reduces to 



Z"(u) + —Z(u) = , 
u 

which has the standard solutions [26] 

Zj{u) = y/uJi(2cjy/u) and Zy(u) = \/uY\(2&-\/u) 
The Bessel functions appearing here have the asymptotic limits 



(3.41) 



(3.42) 



Ji(z) 



cos z 



-7T 



Yi(z) 



2 . / 3 

8111 Z -7T 

7TZ V 4 



which can be used to derive the result 



Z WKB (u) = c+e-^Vd^L [Ji(2w>/u) + iyi(2cjv^)] 



+c_e 



tDu [Ji(2a)\/n) — «Yi(2a)\Af)] 



e " * 



(3.43) 



functions, -ffn 1 ^ = Jn + iY n , H n 2 ^ = J n — iY n , we arrive at the WKB approximation of E 



Recalling finally the relation between Z and E, and using the definition of the Hankel 

— Jn 
,WKB (u) 



a 



1 — u 



(3.44) 



To arrive at an expression for the desired retarded Green's function on the field theory 
side, we still need to expand the WKB solution near u = to obtain an expression similar 
to eq. (3.15), 



E WKB (u) 



—^[l-i—e 2 ^ 110 ) + ulnu(...) 
itujJu \ C-l 



ILO 



IT 



-1 + 27 - m + 21nw) 



2iujao 



(-l + 27 + ivr + 21nw) 



+0(u 2 



(3.45) 
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Figure 6. The behavior of the functions R(uj) (blue curve) and R wkb (oj) (dashed red) of eq. (3.25) 
for r s = 5.01 (smaller amplitude) and 5.5, with = 5. 



From here, we can finally identify A and B and use eq. (3.16) to evaluate 



n WKB M 



N 2 T 2 B N 2 u 2 



8 A 32vr 2 



21nw - 1 + 27 



ITT 



1 + 



1 



A c — p2iu)ao 



(3.46) 



where the ratio c_/c+ is given in eq. (3.12) and the contribution proportional to o) 2 lna) 
can be identified with the T = Green's function [27]. It is remarkable that this result has 
an identical structure with the retarded correlator for a scalar field, cf. eq. (48) of [14]. 

Finally, we note that the large oj limit of eq. (3.19) agrees with the WKB approxima- 
tion, as 

Aout £>^>1 



At 



-ie 



2iuiao 



(3.47) 



an effect demonstrated in fig. 6. From the phase difference of infalling and outgoing waves 
near the horizon one furthermore finds the result 



2iLoao 



exp 



-i& In [ 1 



which allows one to derive an "echo time" in the limit r. 



r h [14] , 



t, 



echo 



1 

2^ 



In 1 



+ 2a 



(3.48) 



(3.49) 



This quantity has a physical interpretation as the time it takes for a wave to travel back 
and forth the WKB potential, and can be identified as the frequency of the oscillations 
witnessed in the function R(u, r s ). 



4 Conclusion 

In the paper at hand, we have used the AdS/CFT correspondence to compute the re- 
tarded Green's function of a U(l) gauge field immersed in strongly coupled, thermalizing 
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N = 4 Super Yang-Mills plasma. On the gravity side, thermalization was modeled via 
the gravitational collapse of an infinitesimally thin spherical shell, which was treated in 
a quasistationary approximation, valid in the last stages of the collapse. The goal of our 
work was to study the production of dileptons in an out-of-equilibrium setting — a quan- 
tity directly proportional to the imaginary part of the retarded correlator as long as the 
Fluctuation Dissipation Theorem is valid. In addition, we studied the flow of the poles of 
the corresponding Green's function, offering an alternative way to follow the thermalization 
process. 

The main results of our work are visible from the figures of section 3.3. The relative 
deviation of the spectral density from its thermal limit, defined in eq. (3.25) and displayed 
in fig. 2, is seen to exhibit fluctuations that increase in frequency and decrease in amplitude 
as the shell approaches the horizon and thermal equilibrium is reached. We also observe 
that the thermalization process is of the top-down type, as the modes with higher energy 
thermalize first. At the same time, figures 3, 4 and 5 display the flow of the poles of the 
retarded Green's function as thermalization proceeds, giving the masses and decay widths 
of the excitations in the system. Three different types of modes were identified from these 
figures, and their physical interpretation discussed. 

All calculations presented in this paper were performed in the quasistatic approxima- 
tion, where the motion of the falling shell was assumed to be slow in comparison with 
the relevant time scales associated with the physical processes studied. To quantify this 
statement, let us note that a characteristic time scale of equilibration can be obtained from 
the functional dependence of the shell radius r s on t. Solving r s from the ds 2 ~ limit of 
eq. (2.1), we obtain an exponentially slow approach towards the horizon [12], 

rs{t) ~ rh = e-'A (4.1) 

where r = l/(47rT). Denoting by oo again the frequency variable in the Green's functions, 
we obtain from here the consistency condition 

<u - ^ « „, (4.2) 



r,(t) 



which is observed to be satisfied for sufficiently high frequencies or late times. Furthermore, 
plugging in eq. (4.1) T > 200 MeV, one obtains the thermalization time scale r < 0.1 fm, 
which one might compare with the typical production time of dileptons with mass/energy 
Q = uj larger than 5 GeV, r p < 0.04 fm. It appears that dilepton pairs produced early on 
have a reasonable chance to escape the system while it is still out of thermal equilibrium. 

While already our present results have provided interesting signatures of dilepton pro- 
duction, we would like to emphasize that this work should be considered only a first step 
towards a more complete holographic treatment of the phenomenon in an out-of-equilibrium 
setting. The next goal would clearly be to try to relax the quasistatic approximation in 
the present setup, taking a step towards fully dynamic thermalization (for closely related 
work, see e.g. [35, 36]). Another assumption made in the current work that would be 
nice to relax is that of working in a spatially homogenous and isotropic background: The 
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holographic system closest to the initial conditions of a heavy ion collision is clearly that 
of two colliding shock waves [37, 38]. In addition to this, it would of course be highly 
interesting to determine other, related quantities, such as prompt photon production, in a 
thermalizing plasma [39]. 
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A Massive scalar field 

In this appendix, we consider the case of a scalar field with mass m, frequency oj and 
vanishing three- momentum in AdS3 space, in order to make sure that our results for this 
somewhat simpler system agree with those of [11, 12]. Introducing the variable u = r 2 /r 2 , 
the equation of motion for the field 4> reads [40, 41] 



lt(l — U)(j) (u>, u) — U(f) (to, u) + 



to 2 rh 2 



1-u 



(co,u) = 0, (A.l) 



where w = oo/{2r} l ) and rh = m /A. Inside the shell, r < r s , the frequency reads again 



tjjin = oj/\/1 — r^/r 2 , and the solution is given by a linear combination of Bessel functions 

J±u [11], 



1 



^inside^, r) = - J u ( — ) - e**» J- v ( ^ ) , (A.2 



T* \ / V IT 

with v = vl + m 2 . 

For the exterior solution at r > r s , we again take a linear combination of infalling and 
outgoing waves, 

^outside (w,r) = c + 4> in (uj,r) + c_0out(w,r) , (A.3) 

where we have defined 

cj) in (u>,r) = z ±l,z '(1 - zf~ 2 Fi(±iiu + f3-;±iuj + ±2i6j; z) , (A.4) 

out 

and (3± = A±/2 with A± = 1 ± \/T+ rn 2 = 1 i v. The hypergeometric function iF\ ap- 
pearing in this expression is regular in the limit r — > r^, in which we obtain the asymptotic 
behavior 

2_ 2 
-r- iu ] n ( r _ r h \ 

cj)- m (w,r -)■ r h ) -> e JaJt e r ft . (A.5) 

out 
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Figure 7. The behavior of R(ui,r s ) for r/j = 5. In order of increasing amplitude, the curves 
correspond to r s = 5.001, 5.01, 5.05, 5.1. 

The coefficients in eq. (A. 3) are finally fixed by the matching conditions of eq. (2.13), and 
hence in the full retarded Green's function on the boundary (r —> oo) 

G R (u, r s ) = G therma V) x H(co, r s ) , (A.6) 



we have 



1 I LUU r> 

G therma V) - " 



H(oj,r s ) 





-to + 


r 2 (- 


-to + 




Sout 

- B in 




>4_out 

' s 



A 



2 



I» 



B rt T(l + 2^^)^ 2 (-to + A 



2 



& ™ r(i - 2iw)r 2 (to + A 



2 



. . r(i + 2to)r 2 -to + 



A ^ r(i - 2to)r 2 (to -j- - 



(A.7) 

In fig. 7, we display the behavior of the relative deviation of the spectral density, R(co, r s ), 
for various values of r s . A behavior qualitatively very similar to the vector case is observed, 
with damped oscillations in the spectral function as well as a top-down type thermalization 
pattern. A study of the flow of the corresponding quasinormal modes (not displayed here) 
is furthermore seen to fully agree with the findings of [11, 12]. 

B Relating the Wightman function to the retarded correlator 

An integral step in our evaluation of the dilepton production rate was to relate the vector 
field Wightman function to the corresponding retarded correlator in eq. (3.2). The Wight- 
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man function can be expressed in terms of the G12 correlator of the real time formalism 
through G < (oo) = ie~ u) /( 27 " , Gi2(oo), 4: and a special case of the Fluctation Dissipation The- 
orem is then required to relate the latter to the retarded Green's function. Though this 
is by far not a trivial identity (for a counterexample, see ref. [36]), we will argue that it 
holds in our quasistatic approximation, where at each stage of its motion the falling shell 
is treated as a stationary object residing at some radius r s > r^. 

To inspect the above statements explicitly, we follow the treatment of a scalar field by 
Herzog and Son [25], modifying it slightly to accommodate the fact that unlike in thermal 
equilibrium, our field exhibits both infalling and outgoing (with respect to the horizon at 
r = r/j) components outside the shell radius. To this end, we generalize eq. (3.12) of [25] 
to read 

Wside(w,r) = ^2 |afc(c+^2,fc + c*_U3 }k ) + /3fc(c + M 4ife + c_-ui ifc ) j , (B.l) 



where we have taken into account that the mode functions 112^ and 113 ^ (for definitions, 
see [25] 5 ) are defined with negative frequency and have subsequently used our eq. (3.13). 
The computation then proceeds in perfect analogy with [25], except that in their eq. (3.16) 
the functions fk are to be replaced by 

4Mi±£y*M. (B . 2) 

c* + + c*_ 

In particular, eqs. (3.17) and (3.18) of [25] are seen to hold, and hence we arrive at their 
relations (2.11)-(2.14) for the Green's functions, where the retarded correlator can be 
evaluated as in our eq. (3.16). 
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